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Abstract 

Most of the empirical studies on stochastic volatility dynamics favor the 3/2 specification 
over the square-root (CIR) process in the Heston model. In the context of option pricing, 
the 3/2 stochastic volatility model is reported to be able to capture the volatility skew 
evolution better than the Heston model. In this article, we make a thorough investigation 
on the analytic tractability of the 3/2 stochastic volatility model by proposing a closed- 
form formula for the partial transform of the triple joint transition density (A, I, V) 
which stand for the log asset price, the quadratic variation (continuous realized variance) 
and the instantaneous variance, respectively. Two distinct formulations are provided 
for deriving the main result. The closed-form partial transform enables us to deduce a 
variety of marginal partial transforms and characteristic functions and plays a crucial role 
in pricing discretely sampled variance derivatives and exotic options that depend on both 
the asset price and quadratic variation. Various applications and numerical examples on 
pricing exotic derivatives with discrete monitoring feature are given to demonstrate the 
versatility of the partial transform under the 3/2 model. 
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1 Introduction 


Stochastic volatility models (SVMs) were introduced to option pricing theory to resolve the 
incapability of the Black-Scholes framework in capturing the volatility smile/skew. Despite 
the difference in the specihc assumption on the instantaneous volatility dynamics, the com¬ 
mon practice of randomizing the volatility is to model the instantaneous volatility/variance as 
a correlated diffusion process. As pointed out by Itkin (2013), classic SVMs use a constant 
elasticity of variance (CEV) process for the instantaneous variance and there are just a few 
choices for the CEV parameter 7 that exhibit mathematical tractability. Nevertheless, various 
versions of SVMs have been proposed in the literature. Hull and While (1987) model the 
instantaneous variance process as a geometric Brownian motion (7 = 1). Scott (1987) and 
Chesney and Scott (1989) let the log variance process be a mean-reverting Ornstein-Uhlenbeck 
(OU) process. Stein and Stein (1991) and Schobel and Zhu (1999) assume that the instanta¬ 
neous volatility follows a mean-reverting OU process (7 = 0). Heston (1993) instead proposes 
a mean-reverting square root process (7 = 1/2) for the instantaneous variance. Barndorff- 
Nielsen and Shephard (2001) use a mean-reverting OU process to model the instantaneous 
variance process. Other variations of SVMs can be found in Bates (1996), Carr et al. (2003), 
Carr and Wu (2003), etc. Among all the proposed SVMs, the Heston model and its variants 
are the most popular ones. It is known that the Heston-type SVMs are essentially a subclass 
of the affine model family [see Duffie et al. (2000)]. 

Despite its popularity, there is a surprisingly large amount of empirical studies that report 
inconsistency of the affine models with market observations. Poteshman (1998) studies S&P500 
index option prices over a 7-year period and concludes that both the physical and risk-neutral 
drifts of the instantaneous variance are nonaffine. Also, the volatility of variance is observed 
to be an increasing convex function of the instantaneous variance. Using an affine drift CEV 
process for the instantaneous variance to £t the S&P500 daily returns over a 30-year period, 
Ishida and Engle (2002) estimate the CEV parameter 7 = 1.71. In a similar work by Jones 
(2003), 7 is found to be 1.33 for daily S&PlOO returns and implied volatilities over a 14-year 
period. Chacko and Viceira (2003) employ the technique of the generalized method of moments 
(GMM) on a 35-year period of weekly returns and a 71-year period of monthly returns and 
estimate the CEV power to be 1.10 and 1.65, respectively. Javaheri (2004) tests three CEV 
power values: 7 = 0.5,1,1.5 on the time series of S&500 daily returns and finds that 7 = 1.5 
outperforms the other two. Using the same data as Jones (2003), Bakshi et al. (2006) test 
several SVMs on the time series of S&PlOO implied volatilities and find that a linear drift 
model is rejected and the coefficient of the quadratic term is highly signihcant. Moreover, 
their estimate of the CEV power is 1.27. The implication of these empirical findings is that 
one should use a diffusion process with nonaffine drift and CEV power greater than 1. With 
a quadratic drift and CEV power equal to 1.5, the 3/2 stochastic volatility model is obviously 
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more firmly supported by the empirical study than the affine models. 

Interestingly, the 3/2 model is not a brand-new model. In fact, it has already been examined 
by Heston (1997) and Lewis (2000) for its analytic tractability, applied to construct short rate 
models by Ahn and Gao (1999) and used to model credit default intensity by Andreasen (2001). 
However, until very recently the 3/2 model has not been considered as a mainstream SVM. The 
growing attention from the academia to the 3/2 model is partially attributed to the increasing 
interest in consistently modeling equity and volatility markets (represented by VIX), as a result 
of the booming of the volatility derivative market. In a recent work by Carr and Sun (2007), 
a new framework for option pricing that directly models the variance swap rate is proposed. 
They argue that the 3/2 specihcation for the instantaneous variance is a direct consequence 
of the model consistency requirement. Following Carr and Snn’s work, Itkin and Carr (2010) 
introduce the 3/2 power time change process, and Chan and Platen (2010) price long dated 
variance swaps nnder the 3/2 stochastic volatility model. By performing extensive nnmerical 
comparisons between the Heston model and the 3/2 model, Drimus (2012) reports that the 3/2 
model is snperior to the Heston model in the sense that it is able to predict upward-sloping 
volatility of variance smiles, which is in good consistency with market observations. Most 
recently, Goard and Mazur (2013) report strong empirical evidence that VIX follows a 3/2 
process other than an affine square root process. In an effort to consistently modeling VIX and 
eqnity derivatives, Baldeaux and Badran (2014) consider a 3/2 plus jumps model for pricing 
VIX derivatives. Ynen et al. (2014) derive closed form pricing formnlas for exotic variance 
swaps nnder the 3/2 model. 

Analytic tractability of an nnderlying model is essential for derivatives pricing. As for 
the 3/2 model, the characteristic fnnction of the log asset price process is derived in Heston 
(1997) and Lewis (2000). In order to facilitate the pricing of volatility derivatives, Carr and 
Sun (2007) obtain the joint characteristic fnnction of the log asset price and the quadratic 
variation, and Baldeaus and Badran (2012) extend their result to the 3/2 plus jumps model. 
Lewis (2014) derives the joint transition density of the log asset price and the instantaneous 
variance for the 3/2 model with constant parameters. Nevertheless, the full characterization 
of the 3/2 model seems to rely on the joint distribntion of the triple {X,I,V) which stand 
for the log asset price, quadratic variation and the instantaneous variance. Our paper hlls 
this gap by providing a complete description of the joint distribution through the closed-form 
partial transform of the triple transition density. Two distinct formulations are proposed for 
deriving our main result. The partial differential eqnation (PDF) approach is inspired by the 
work by Carr and Sun (2007) and Lewis (2014). Carr and Snn (2007) solve the governing PDF 
with an ingenious choice of substituting variable and then reduce the PDF to an ODF. Lewis 
(2014) uses a sequence of transformations to reduce the original governing PDF to a first-order 
linear PDF which is solvable by the standard method of characteristics. Motivated by the joint 
exponential affine structure of the log asset price and qnadratic variation, we dehne the partial 
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transform of the triple joint density function and then solve the governing PDE by converting it 
to a Riccati system of ordinary differential equations through a Laplace transform. In the pure 
probabilistic formulation, the partial transform is factorized as a product of the conditional 
characteristic function of the integrated variance and the marginal transition density of the 
instantaneous variance. We then End the explicit expressions for both terms by using the 
change of measure technique and the reciprocal relation between the 3/2 process and CIR 
process. 

Put in the relatively thin collection of literature on option pricing under the 3/2 model, 
the contribution of our work is three-fold. First, we propose a closed-form formula of the 
partial transform of the triple transition density under the 3/2 model with a time-dependent 
mean reversion parameter. The closed-form partial transform fully characterizes the analytic 
property of the 3/2 model, and most of the existing analytic formulas for the joint density 
functions and characteristic functions in the literature can be viewed as marginal versions of the 
newly derived partial transform. The closed-form partial transform enables us to accommodate 
the standard transform method to the pricing of derivative products whose terminal payoffs 
have exotic dependency on the asset price and quadratic variation under the 3/2 model with a 
time-dependent mean reversion parameter. Second, we provide two distinct derivations of the 
main result and each approach exhibits some unique novelty and improvement over the extant 
methodologies. In particular, our PDE approach takes advantage of the affine property of the 
pair {X, I) and explicitly solves the governing PDE by converting the equation into a Riccati 
system of ODEs. Our probabilistic approach explores the relationship between the partial 
transform and the conditional characteristic function of the integrated variance and works out 
the closed-form expression for the latter with the change of measure technique. Third, using 
the newly derived formula, we manage to analytically price a variety of lately developed exotic 
derivative products, including finite-maturity discrete timer options and discretely sampled 
weighted moment swaps. 

The rest of the paper is structured as follows. In Section 2, we provide a detailed description 
of the 3/2 model and its financial intuition. Section 3 is devoted to the derivation of our 
main result on the partial transform of the triple joint transition density under the PDE and 
probabilistic formulations. We also show that most marginal characteristic functions under the 
3/2 model can be derived immediately from the partial transform of the triple. Section 4 covers 
examples for demonstrating the applications in derivatives pricing under the 3/2 model. In 
particular, the pricing of finite-maturity discrete timer options and discretely sampled weighted 
moment swaps are investigated in details. Numerical experiments and analyses are provided 
in Section 5 and conclusive remarks are given in Section 6. 
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2 The 3/2 model 

Consider the 3/2 stochastic volatility model specihed as follows: 

^ = (r - 5) d« + + yrVdM-?), 

( 2 . 1 ) 

dCt = Vt{et - KVt) dt + ml 

where hh/ and are two independent Brownian motions nnder the pricing measnre Q. Here, 
we assnme a constant riskfree rate r and dividend yield q. Time-dependent deterministic 
riskfree rate and dividend yield can be accommodated withont difficulty. In contrast to the 
square root process with affine drift in the Heston model, the parameters in the 3/2 variance 
dynamics need to be interpreted differently. The speed of mean reversion is now kV), which 
is linear in V) and is a stochastic quantity. Since k > 0 under usual scenarios, the mean 
correction is quicker when the instantaneous variance is higher. Also, £ cannot be interpreted 
as the same volatility of variance as in the Heston model. In fact, one needs to multiply it by 
a scaling factor V) in order to make it comparable to its counterpart in the Heston model. The 
long-term mean reversion level is given by Ot/n. As pointed out by Itkin and Carr (2010), 9t 
could be an independent stochastic process in the most general setting. By conditioning on 
the path of 0^, the analytic tractability of the 3/2 model with stochastic mean reversion level 
remains intact. Although in practice 9t is often taken to be constant, this flexibility allows for 
more delicate modeling of the instantaneous variance dynamics when necessary. Thus, in this 
paper, we assume 9^ to be a deterministic function of time, but an extension to an independent 
random process 9t is not difficult. 

Similar to the Heston model, not all choices of model parameters are admissible, in the 
sense that the non-explosion of V) and the martingale property of the discounted asset price 
process are guaranteed. According to Drimus (2012), the parameters of the 3/2 model specified 
by 112.11) are constrained by 

K-pe>-^. ( 2 . 2 ) 

Notice that under normal market condition, k > 0 and p < 0, the above inequality is auto¬ 
matically satisfied. 


3 The closed-form partial transform 

In this section, we present our main result on the closed-form formula for the partial trans¬ 
form defined previously. On top of deriving the main result by solving the governing PDE 
analytically, we provide an insightful probabilistic argument and some interesting intermediate 
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results are obtained as a byproduct. We also briefly discuss how to deduce various marginal 
characteristic functions and transforms from the main result. 


3.1 The main result 


Let Xt = In St be the log asset price and It = Vg ds be the quadratic variation of the log asset 
price process. Let G{t,x,y,v;t',x',y',v') be the joint transition density of the triple {X,I,V) 
from state {x, y, v) at time t to state (x', y', v') at a subsequent time t'. By the Feynman-Kac 
Theorem, G satisfies the following Kolmogorov backward equation: 


dG 



^ v d^G 

dx 2 dx"^ 


dG ^dG 

+ V— + v[0t - K,v)— + 
dy dv 


di^G 
2 dv"^ 


+ pev^ 


d^G 
dxdv ’ 


(3.1) 


subject to the terminal condition: 


G{t\ X, y, V, t', x', y', v') = d(x - x')d(y 


y')6(v-v'), 


where 6(-) is the Dirac delta function. 

The partial transform of G, denoted by G, is defined by 


pco poo 


G{t,x,y,v,t\uj,r],v’)= / x, y, u; F,x', y', u') d|/'dx'. (3.2) 

J —oo J 0 


The reason for introducing G is that the affine structure of the pair {X, I) makes G more 
tractable than G itself. Also, G is more convenient to be used in option pricing. Obviously, G 
also satisfies (I3.1|l with the terminal condition being G{t', x, y, u; f', cu, p, v') = — v'). 

Since (13.ip has no coefficient involving x or y, it is natural to come up with the following 
solution form: 


G{t,x,y,v;t',u,p,v') 




u,r), v'), 


(3.3) 


where g satisfies the following PDE: 


dt 


“(’'■«■ i) 



+ ipv 


g + [v{et - Kv) + iupev'^]^ + 


g2^3 Q2g 

2 dv"^ ’ 


(3.4) 


subject to the terminal condition: 


g{t',v;t',uj,p,v') = 6{v - u')- 


From fl3.4p . we see that g is essentially a function of v and t, while all the other variables affect 
g either through the PDE coefficients or boundary condition. 


Theorem 1. Under the 3/2 stochastic volatility model specified by the partial transform 

of the triple joint transition density function defined by O.^l) is found to be 


G{t,x,y,v,tjuj,p,v’) 


^iujx+iT]y 


g{t,v;tju,p,v'), 


(3.5) 
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where 


g{t,v;t\u,r],v’] 


= e“(*'-*)^exp 
Ct 


+ 1 A /a 

Ctvv' j iv'Y V \Ct\ vv' 


(3.6) 


Here, 


a = iuj(r — q), k = k — iupe, 

, , 1 K V iu + — 2m 

c=Wl5 + jjW 


Proof. We prove it by explicitly solving (13.dh . The details are given in Appendix El □ 


Remark 1. The solution form of G given by (13. 3 p holds for any stochastic volatility process. 
In fact, if we consider a general stochastic volatility model specified by 


-J. = (r-q)dt+ ^/V,{pdWl + Pl-p‘dWP, 

Jt 

dV,=a{t,V,)dt + l3{t,V,)dW;. 


(3.7) 


It can be shown in a similar manner that the governing equation for g is given by 


dt 


.“(’'■«■ i) 



+ igv 


.... .,dg /3{t,vYd‘^g 

g + la(t,v) + iwpv^/3(«,ii)]^ + —-— 


(3.8) 


subject to the terminal condition: g{t',v]t',u,g,v') = 6{v — v'). The partial transform of the 
model (13.7p admits a closed-form expression if (13.8p can be solved explicitly. 


Remark 2. Theorem Q can be easily extended to incorporate jumps in the asset price process. 
The analytic formula under the 3/2 plus jumps model is given in Appendix\M 


A probabilistic formulation of the partial transform 


The partial transform dehned by fl3.2p admits an interesting conditional factorization, from 
which we can derive our main result in a pure probabilistic manner. Let py(V)/|ld) be the 
transition density function of the instantaneous variance for any t' > t under the pricing 
measure Q. Since G{t, X*, A, Vp t', Xt>, It>, W) = G{t, Xt, A, Vp t', Xt>, It'\Vt')pv{Vt'\Vt), we have 


G{t,Xt,It,Vt;t',uj,g,Vt/ 


W, A, Id; t', x', y', W) dy'dx' 




(3.9) 
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where Ej[-] is short for E[-|J^t] under the measure Q with {J^t)t>o being the natural hltration. 
The dynamics in fl2.ip can be rewritten to give [e.g. Zeng et al. (2013)] 


Xt' — Xt + (r — q){t' ~ ~ 


In Vt> - In Vt 


6. ds 


+ 


pe 


K 1 
“W + 


Kds+ 


_ J t 


With this expression, the conditional expectation can be reformulated as 





E( 

'^ico{Xt,-Xt) 

[ Vds,Ve 


Vr 




Jt ^ 




_ ^iuiXt+irilt^iui(r-q){t'-t) 





V,> 


.(3.10) 


Therefore, the original problem boils down to the computation of the transition density of 
the variance process PviVt'lVt) and the conditional characteristic function of the integrated 
variance in the form of E[e'^'^‘ Vt]. 

Transition density of the instantaneous variance 

It is well known that the reciprocal of the 3/2 process is a CIR process. Write Ut = then 
Ut is given by the following inhomogeneous CIR process: 


dU, = ((k + £=) _ it - es/UidW^. 


(3.11) 


Lemma 1. For any t' > t, Uf follows a non-central chi-square distribution conditional on Ut- 
Its transition (conditional) density function is given by 


PuiUt'\Ut) = ^exp 


AtUt' + Ut 

Ct 


AtUt> 

Ut 


1 

2 






(3.12) 


where At and Ct are given in TheoremUl 


Proof. The inhomogeneous CIR process fl3.11l) equals a squared Bessel process B(t) trans¬ 
formed by the following space-time changes in distribution 


{Ut,t > 0) 


1 

e/o 


B 






Based on the above relation and transition density of the squared Bessel process, Jeanblanc 
et al. (2009) provide the transition density of the inhomogeneous CIR process. We refer 
interested readers to Jeanblanc et al. (2009) for more details. □ 























It follows that the transition density of Vt can be inferred to be 


Pv{Vf\Vt) = ^exp 


AtVt + VA 1 fAtVty^^ (2 

CtVtV, ) V Vf ) \cyvt,Vt) • 


(3.13) 


Conditional characteristic function of the integrated variance 

Because of the strong path dependency of V)], the evaluation of this conditional 

expectation poses a major challenge. Fortunately, with the change of probability measure 
technique, we manage to hnd the analytic formula for this conditional expectation. 


Theorem 2. Under the 3/2 stochastic volatility model specified by H2.1\) . the characteristic 
function of the integrated variance conditional on F)' admits 


E 






2k 

7 ^ 


Ct 





where the notations At and Ct are the same as those given in TheoremUl 


(3.14) 


Proof. The proof relies on a smart choice of a new probability measure under which the evalu¬ 
ation of expectations can be done conveniently. The details can be found in Appendix O D 

Remark 3. Formula fl3.14p can be viewed as a generalization of Baldeaux (2012) for a time- 
dependent mean reversion parameter Ot- Baldeaux’s derivation is based on Lie symmetry and 
Laplace method, whereas in this article we manage to propose an easier alternative to obtain 
a more general result. Furthermore, one can now construct the exact simulation of the asset 
price process under the 3/2 model with a time-dependent mean reversion parameter. 

One can now easily verify that combining fl3.9p . (13.1 Op . fl3.13p and fl3.14p immediately leads 
to fra . 


3.2 Marginal characteristic functions and transforms 

From the partial transform of the joint transition density of the triple {X, I, V), we can deduce 
many useful marginal characteristic functions and transforms, including the renowned formulas 
proposed by Carr and Sun (2007) and Lewis (2014). 

Corollary 1 (Carr and Sun, 2007). Under the 3/2 stochastic volatility model specified by 112.1\) . 
the joint characteristic function of the pair {X, I) conditional on Xt is given by 

Vp, t’, U, T]), (3.15) 
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where 

= M (aj,-^y (3.16) 

Here, 

a = —- —~— + c, (3 = l + 2c, 

I e‘^ 

r is the gamma function, and M is the confluent hypergeometric function of the first kind. 

Proof. See Appendix [Dl □ 

Remark 4. The marginal characteristic functions of X and I can be deduced by simply setting 
Tj = Q and oj = Q, respectively. 


Denote the joint transition density fnnction of {X,V) by Gi(t,x,v,t',x',v') and define its 
partial transform by 


Gi{t, x,v;t',u,v') = I Gi{t,x,v;t', x',v') dx' 


(3.17) 


Corollary 2. Under the 3/2 stochastic volatility model specified by H2.1\) . the partial transform 
defined by ^3.17^ is given by 


Gi{t,x,v]t',u,v') = e^‘^'^gi{t,v]t',u,v'), 


(3.18) 


where 


gi{t, V] t', u, v') = g{t, V] t', u, 0, v') 

, At ( AfV + v' 




Gi 


exp 


Gtvv' 


. l \-2 


Atv 


'2c 



with 



Proof. By virtne of fl3.2p . Gi can be obtained by simply setting 77 = 0 in fl3.5|) . □ 

Remark 5. One can immediately see that this result can be viewed as a generalization of Lewis 
( 2014 ) for a time-dependent mean reversion parameter 9^. 

Remark 6. By further setting cu = 0 m we can retrieve the transition density ofV: 

Gv{t, V] t', v') = g{t, v] t', 0 , 0 , v'). (3.19) 

This hence provides an alternative derivation of LS.13\) . 
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The joint characteristic function of {X, I) given by fl3.15p is useful for pricing European 
contingent claims that only depend on the values of {Xt, It) at the maturity T. However, for 
more exotic derivatives, fld.lSp alone is not enough. For instance, the pricing of discretely mon¬ 
itored hnite-maturity timer options also requires the joint characteristic function of {Xt.,!^.), 
where i ^ j- The evaluation of discretely sampled weight moment swaps needs the joint char¬ 
acteristic function of {Xt-, Xtj, Xt^.), where i,j and k are all different. These multivariate joint 
characteristic functions cannot be deduced from Corollary [T] or [21 and one needs to make use 
of the full characterization of the triple. 

Let uj = (cji, • • • and rj = {rji,--- , 1 ]^)'^ be scalar vectors in M™, and X 7 - = 

(Xt-^, ■ ■ ■ ,Xt^)'^ and I 7 - = (4^, • • ■ be vectors of log asset price and quadratic varia¬ 

tion, respectively, observed at a sequence of ascending future instants denoted by T = {tiljhi- 
Then, we define the multivariate joint characteristic function by 

<!>{t,Xt,It,Vt-,r,u:,ri) = (3.20) 

where t denotes the current time and t <U for any i. Specihcally, the bivariate joint character¬ 
istic function with m = 2 is of particular interest and will be used repeatedly in the subsequent 
section. 

Corollary 3. Under the 3/2 stochastic volatility model specified by H2.1\) . the bivariate joint 
characteristic function defined by l\3.2(j\) with m = 2 is given by 

^{t,Xt,It,Vt,T,oj,'n) 

poo ^3 21 ) 

^ ^i{ui^+uj 2 )Xt+i{m+r, 2 )It / g{t^Vt,tl,OJi+ U 2 ,rii + ri 2 ,v')h{ti,v'-,t 2 , 0 J 2 ,ri 2 )dv', 

Jo 

where g and h are given by H3.6\) and H3.16\) . respectively. 


Proof. By direct calculation, we have 

(^{t,Xt,It,Vt-,T,uj,r}) 

= Corollary Cl) 

poo 

= / h{ti,v'-,t 2 ,U 2 ,r] 2 ) 


f*oo poo 


' —00 Jo 


/U^+^2)x'+i{m+V2)y' Xt, It, Vt; p, x', y', v') dx’dy’dv’ 
= I h{ti,v'-,t2,U2,r]2)G{t,Xt,It,Vt,ti,Ui+U2,r]i + r]2,v')dv' 

Jo 

poo 

^ y(0Ji+uJ2)Xt+i{m+V2)It / g{t^Vt;ti,Ui+UJ2,T]l+T]2,v/h{ti,v';t2,UJ2,T]2)dv'. 


This completes the proof. 


□ 


For notational convenience, in the sequel, we suppress the dependence of $ on (Xt, fi, Vt) 
and write d>t{T, UJ, rf) for short when there is no ambiguity. 
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4 Applications in derivatives pricing 


In this section, we demonstrate how the partial transform of the triple and its induced marginal 
characteristic functions and transforms can be used to price various exotic derivative products 
that may be embedded with path-dependent feature, and have a complicated terminal payoff 
structure that depends on the asset price and quadratic variation. 


4.1 Finite-maturity discrete timer options 


The hrst presence of timer option in the literature dates back to the 90s when Neuberger (1990) 
and Bick (1995) discuss the pricing and hedging of such an imaginary product. Since the official 
launch of timer option by Societe Generale in 2008 [see Sawyer (2007)], research work that 
examines this exotic product from either an analytic or approximation perspective has been 
extensive [see Zeng et al. (2013) and references therein]. Voluminous as the literature is, most of 
the existing results are achieved under the assumption of perpetuity and continuous monitoring. 
Direct investigation on the discrete timer option with hnite maturity is relatively rare and is 
more mathematically challenging. In this subsection, we will discuss how to analytically price 
discrete timer options with hnite maturity under the 3/2 model based on our newly derived 
partial transform of the triple density. 

Instead of having a deterministic expiry date as in a vanilla option, a discrete timer option 
with hnite maturity expires on a random date which is either the hrst time when a pre-specihed 
variance budget is fully consumed by the realized variance of the underlying asset price or the 
pre-specihed mandatory expiry date, depending on which one comes earlier. Let T be the 
mandatory expiry of the timer option and denote the tenor of the monitoring dates for the 
realized variance by {foGi ‘ ‘ ‘ Ga^}- For brevity, we assume equally spaced monitoring interval, 
i.e. tj = jA = for j = 1, 2, • • • , iV. 

At the initiation of the timer option, the investor specihes an expected investment horizon 
To and a target volatility uo to dehne a variance budget 

B = a^To. 


Let tb be the hrst time in the tenor of monitoring dates at which the discrete realized variance 
exceeds the variance budget B, namely. 


Tb 


mm < j 




A. 


Let the current time be to = 0, then the price of a hnite-maturity discrete timer call option can 
be decomposed into two components, depending on whether the variance budget is consumed 


12 



up before the mandatory maturity or not. 


Co = Eo[e-'-(^^"«)max(^TAr,-i^,0)] 

= Eo[e"’’'^max(S'T - K, 0)1{^^>t} + e"''^®max(S'^^ - K, 0)1{^3<t}], (4.1) 


where K is the strike price. 


To explicitly calculate the timer option price, we use the quadratic variation It as a proxy 
of the discrete realized variance for the monitoring of the hrst hitting time. That is, we redehne 
tb as follows: 


Tb = mm < j 


4 >B'>A. 


With this simplihcation, the price of a hnite-maturity discrete timer call option can be conve¬ 
niently computed by further decomposing it into a sequence of timerlets as follows: 


Co = E[e ''^max(S'r - K,0)1{i^^b}] 


E 


N-l 


e - K, 0)1 {b.<b} - max(4,^^ - K, 


L j=o 


(4.2) 


The key observation is that the event {tb > t} is equivalent to {It < B}. That explains the 
equivalence between the hrst term of fl4.1l) and fl4.2p . Other terms can be deduced similarly 
by noting that 


N-l N-l 

{tb < r} = U [U <tb < U+i} = 

i=0 i=0 

The above decomposition is hrst proposed by Lee (2013) and has been applied in Cui (2014). 

The series of expectations in fl4.2p can be easily evaluated by the standard transform method 
[e.g. Lee (2013)], if we know the explicit expressions of the joint characteristic functions of 
(Wj, 4) and (Xt^,^j,4)- A direct application of Corollary [1] and Corollary [3] gives 


poo 

Eo[e“^4+i+'^4] = / g(^Q^ Y^. v')h{tj, u'; fj+i, w, 0) du', 

Jo 

Note that the Fourier transform of ( 4+1 — and ( 4+1 ~ admit the 

same analytic representation 


F(u;,p) 






K)^l{y<B} dxdy 


(iw -f- a;^)4’ 
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where uj < —1 and rjj > 0. By Parseval’s theorem, the hnite-maturity discrete timer option 
price can be derived as follows 


^OO POO 


Cn = 




dTT^ 


ii^Xt +ir)7t 


' —CXD J —OO 


N-l ^ 
j=o 


OO POO 


e-rtj+1 


' —OO J —OO 


-| ^OO POO 

/ F[b:,r])H{b:,r])dib:nriR, 




' —OO O —OO 


where 


N-\ 


H{u, 7]) = H; tv, w, v) + ^ 

j=0 

f *00 

g{0, Vo] tj,u, 7], v')h{tj, v'] tj+i,u, 0) dn' - h(0, Vq] tj+i, u, r]) 


4.2 Discretely sampled weighted moment swaps 

Let 0 = to < ti < • • • < tv = P be iV + 1 sampling dates and {4}^^ be a sequence of integers 
that take value from the index set {0,1, 2, • • • , N}. We dehne the kth weight Wk as a mapping: 
Sti^ I—t for some function /. Then, the floating leg of a weighted mth moment swap is 

dehned by 

Immediately, we identify the following special cases of discretely sampled weighted moment 
swaps. 


product 

ik 


fix) 


m 

variance swap 

- 


X 

1-4 

1 

2 

gamma swap 

'ik ~ 

= k 

X 

I-)- 

x/So 

2 

corridor variance swap 

ik ~ 

= k OT ik = k — 1 

X 

H- 

^{l<x<u} 

2 

self-quantoed swap 

'ik ~ 

= N 

X 

1-4 

x/So 

2 

skewness swap 

- 


X 

1-4 

1 

3 


Table 1: Various weighted moment swaps 


The pricing of the above weighted moment swap requires the computation of the fair strike 
price K which is set to be 




■ V 

E 

.fc=i 




In 








^k — l 
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such that it costs zero for both parties to enter into the swap contract. Due to additivity of 
expectation, we have K = ^ 'l2k=i where 


Lh — En 


fiSuJ In 


5. 






^k — 1 


To compute Lk, let us hrst consider the trivial case where ik = 0 or f{x) = 1, meaning that 
the weight is deterministic. As a result, we have 


Lk — WfcEo 

= WkEo 
Wk & 


In 




St,., 

1 a™ 


jm 




0=0 


im 

Wk r 

im 




b=0 


Gy (to, Do; 4-1, n')h(4_i, 4; tk, </>, 0) dn' 


b=0 


The above calculation applies to variance swap and skewness swap. Note that in the last 
but second equality, the expectation is identihed as the forward characteristic function hrst 
proposed by Hong (2004) and later used by Itkin and Carr (2010) to price discrete variance 
swaps. 

Now, suppose ik ^ 0 and / is also nontrivial and assume that / admits the generalized 
Fourier transform with respect to In as follows 

/ CXD 

/(e")e-""da;, 

•OO 


where the transform variable a; is a complex number and its imaginary part uj is hxed in a 
way such that |/(ci;)| < oo. Then, we have 


Lk — En 


1 . 1 

— / duR - piHXt,-Xt^_d 

271 


27ri’^ 




jm Q(j)m 


b =0 


0=0 
dcui?. 


(4.3) 


We remark that in fl4.3p . Eq 
terms (including the differen 


+i 0 (Xtj, ^ model-specihc term, whereas all other 

dal operator) can be considered as product-specihc. This is a 
standard integral representation with the Fourier transform method, which holds for any pricing 
model. 


In order to calculate the expectation in fl4.3l) under the 3/2 model, we hrst consider the 
scenario where 4 ^ k. Using the tower rule and the bivariate joint characteristic function 
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fl3.2ip . we obtain 


Eo 
= Eq 

= Eq 






gi(a;+fli)Xt^_^ O)5(i(4_i,i/4_^;4,a; + 0,n')dn' 


(*CO /•CO 


_ giojXo 


'0 ^0 


h(tk, v'; tii^,u, 0)gi{tk-i, v; 4, uj + (f>, Vq; 4-i, tu, v) dn'dn. (4.4) 


Specifically, when 4 = k, the double integral can be further simplihed as a single integral. 
Since h(tk, v'; tk,u;, 0) = 1 and h(tk-i, v; tk, oj + (j),0) = gi(tk-i, v; tk, w + 0, v') dn', we have 


En 


^iujXti^ +i4>(Xti^ ) 


_ ^icdXo 


h( 4 _i, v; tk, uj + (j), 0 ) 5(1 (4, Vq; tk-i,uj, v) dn. 


In a similar manner, we obtain the following formula for the case 4 4^ — 1: 


En 


^iojXt +i0(Xtj, - Vtj._ J 
^co poo 


— gii^^o 


poo 

h(tk-i,v'-, tk, 0, 0)Gv(tii^,v; 4-1, v')gi(to, Vq; ti^,uj, v) dn'dn. (4.5) 


'0 Jo 


Further simplihcation is possible when ik = k — 1. In that case, Gv(tk-i,v,tk-i,v') = 6{v — v'), 
and hence we have 


En 


+i(l>{Xtf. - Vt ) 


_ gitjXo 


h(4_i, v; tk, 0, 0 ) 51 ( 4 ,14; 4 - 1 , 1 :^, v) dn. 


Observe that both fl4.4p and fl4.5l) indicate that only one function in the integrand depends 
on the dummy variable 0. As a result, the differentiation with respect to 0 in fl4.3p can be 
conveniently performed on that single function. The hnal formula for Lk is obtained as a triple 
integral (a double integral in the two special cases) by plugging fl4.4|) or fl4.5l) into fl4.3p . 

Self-quantoed variance swaps 

As an illustrative example, we consider the fair strike price of a self-quantoed variance swap. 
According to Table (H its fair strike is computed by 



Apparently, Lk is now determined by fl4.4p . Moreover, 


/(^) 


1 



gXg-ia;x 


dx 


2x5{uj + i) 
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Consequently, the integration in fl4.3p can be computed explicitly as follows: 


Lh 


1 


En 




6=0 


poo /*oo 


Jo Jo 
Therefore, we have 


h(4, v'; T, -i, 0)—5fi(4-i, v; 4, 0 - i, v') 4-i, -i, v) dn'dn. 


0=0 


POO POO ^ 

K = -- / ^h(4,4;T,-i,0) 

J- Jo Jo 


90^ 


(4.6) 


Sl(4-i> f; 4, *6-6 1*') 51 (4> Ki; 4_1, -i, u) du'dn. 


0=0 


5 Numerical examples 

In this section, we show some numerical examples that illustrate the performance of pricing 
formulas of the timer option and variance derivatives under the 3/2 stochastic volatility model. 
We also examine the pricing behaviors with respect to varying model parameters. 

Table [2] lists the default parameter values of the 3/2 model in our sample calculation. These 
parameter values, which have been calibrated to the S&P 500 option prices, are taken from 
Drimus (2012). All our numerical examples will use this set of parameter values unless specihed 
differently. 


K 9 

£ W 

P 

5o 

r 

q 

22.84 4.979 

8.56 0.060025 

-0.99 

100 

0.015 

0 


Table 2: Parameter values in the 3/2 stochastic volatility model 

Finite-maturity discrete timer options 

We set the strike price to be iC = 100, and the number of monitoring instants to be W = 100. 
We also reset the value of correlation coefficient p = —0.5 and initial variance Vq = 0.087 which 
are adopted in Zeng et al. (2013) for comparison purposes. Considering that extensive analysis 
has already be done in Zeng et al. (2013) and references therein, here we mainly explore the 
impact of variance budget on finite-maturity discrete timer option call prices over a wide range 
of maturities. 

Figured] shows that the finite-maturity discrete timer option price is an increasing function 
of both the maturity T and the variance budget B. Intuitively, an option is usually more 
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expensive with a longer exercise time. The higher variance budget B or the longer the maturity 
T leads to the later arrival of the exercise, thus giving a more expensive timer option price. 
When T becomes sufficiently large, the hnite-maturity discrete timer option value becomes 
almost insensitive to T and converges to that of the perpetual discrete timer call option from 
below. The hnite-maturity discrete timer option value exhibits a more pronounced convergence 
given a lower B. 



Mandated maturity 

Figure 1: Plot of the hnite-maturity discrete timer call option price versus mandated maturity 
given diherent values of variance budget. 

Self-quantoed variance swaps 

We perform numerical analysis on the self-quantoed variance swap. Specihcally, the fair strike 
prices of the self-quantoed variance swaps are compared with those of the variance swaps to 
exemplify the unique features of the former. In Figured two plots of the fair strike price versus 
the correlation coefficient p and volatility of variance e are given for half-year daily sampled 
swap contracts {N = 126, At = 1/252). 

In Figure [2] (a), we vary the volatility of variance parameter while hxing other parameters 
at the values given by Table [21 to analyze the sensitivity of the fair strike prices of the self- 
quantoed variance swap and the vanilla variance swap to £. Within a realistic scope of e, we 
observe large deviation in the fair strike prices of both swap products. Moreover, the fair strike 
prices are decreasing in e. An intuitive explanation for this phenomenon can be found in Yuen 
et al. (2014). By contrasting the self-quantoed variance swap to the variance swap, we see that 
the disparity between the two gets wider as £ decreases. That’s because the protective feature of 
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the self-quantoed variance swap begins to reinforce when the realized variance is getting large. 
A similar sensitivity analysis with respect to p is performed in Fignre [2] (b). As expected, 
the self-qnantoed variance swap is mnch more dependent on the correlation coefficient than 
the variance swap. When p < 0, we have the renowned leverage effect and the self-qnantoed 
variance swap is embedded with a crash protection for the seller. On the other hand, when 
p > 0, the soaring realized variance may be farther exacerbated by the rising asset price, 
resnlting in a larger risk exposnre. 


(a) fair strike versus voi of variance (b) fair strike versus correiation 



Fignre 2: Comparison of the fair strike prices of a half-year daily sampled self-qnantoed variance 
swap and a vanilla variance swap, with varying valnes of (a) vol of variance, (b) correlation 
coefficient. 


6 Conclusions 

We explore the analytic tractability of the 3/2 stochastic volatility model by investigating the 
joint transition density of the triple consisting of the log asset price, qnadratic variation and 
instantaneons variance. We obtain the closed-form partial transform of the joint triple density 
nnder the 3/2 stochastic volatility model with a time-dependent mean reversion parameter 
nnder two different formnlations. In one approach, we establish the governing PDF of the 
partial transform which is then transformed to a Riccati system of ODEs and solved explicitly, 
while in the other one we relate the partial transform to the conditional characteristic fnnc- 
tion of the integrated variance and compnte the latter one nsing a set of probabilistic tools, 
snch as the change of measnre techniqne. The newly derived partial transform serves as the 
nnderpinning of the transform methods for derivatives pricing nnder the 3/2 model and most 
extant characteristic fnnctions or transforms can be viewed as marginal versions of the partial 
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transform. As illustrative examples, the pricing formulas for the hnite-maturity discrete timer 
options and self-quantoed variance swaps are expressed in terms of the partial transform and 
its induced marginal joint characteristic functions. 
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Appendix A Proof of Theorem [T] 


It suffices to solve fl3.4p which can be rewritten as 


lU + UJ 


— yq 


vg = 0, 


(A.l) 


where g = ge°‘^^~^'\ The terminal condition remains to be g{t',v;t',u,q,v') = 5{v — v'). To 
solve flA.ip . we first change the variables (n,n') to {u,u') with u = 1/v and u' = 1/v'. We then 
define f(t,u;t',u;,q,u') that satisfies 


u 

g{t, V-1', u, q, v') = U] t', u, q, u'). 


By expressing the partial derivatives of g in terms of the partial derivatives of /, we obtain the 
governing equation for / as follows 


K 

dt 




r ^2 


+ 


[O' T O') T KO! — 


ica + oj^ 


+ iq 


/ 

u 


subject to 

fit', u] t', UJ, q, u') = 5{u — u'). 
Apparently, if we choose a = a{uj,q) such that 


-^yCX -|- (T) "t~ KCX 


iu; + 


+ iq = 0, 


then the governing PDE of / becomes 


df , r 2/ , ~ n n y 

(a + l) + K-«,«]^-ae,/, 

where all the coefficients are affine in u. It follows that a can take two values 

where 


, , , ,T k\^ iuj + uj'^ — 2iq 

c = c{uj,q) = \l(- + ^] + 


(A.2) 


(A.3) 


There are various ways of solving flA.2D . We hnd the following approach that makes use of the 
Riccati system of ODEs the best. An alternative method is appended to the end of this proof. 


Riccati system of ODEs 
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We write the Laplace transform of / with respect to u' as follows 

poo 

j ,u:,7],u')(lu'. (A.4) 

Then, / satishes 

subject to the terminal condition f(t',u;t',u;,ri,^) = e““^. By virtue of the affine structure of 
the PDE coefficients, flA.Sp admits the following exponential affine solution: 


/(f,M;f',a;,? 7,0 = exp{B{t,^)u + D{t,^)), 


(A.6) 


where and D{t,^) are parameter functions determined by the following Riccati system 

of ODEs: 


dt 2 ' ’ 

dD 

= [5: (oTijTKj-B — caOti 


dt 

with boundary conditions B{t',^) = and D{t',^) = 0. It can be found that [e.g. Zheng and 
Kwok (2014)] 

B{t,0= ^ 


where 


By noting that 


we obtain 


At + Ct^ 


-2 rt' 


(A.7) 


^ ./f 


d\niAt + CtO . e 

— = —ot — 


dt 


2 Ai + Ct^ 


a + 1 + -^ 


D{t,0 = -^ 

Plugging flA.7p and flA.Sp into flA.Oj) yields 

f[t,u]t = A^ " exp 

Inverse Laplace transform 


\n.{At + CtC) + 




At + ae 


2k 

a 2 -\—- 


{At + Ct^) 


Og ds. 


-2a-2- 


2 ft 

72 - 


(A.8) 


(A.9) 


Finally, / can be retrieved by taking the inverse Laplace transform of / with respect to f. For 
h{p) = p~^~^e'^/^ for V > —1, the inverse Laplace transform of h{p) is given by 

pT+ioo / \ul 2 

= e^^h{p)dp=^-j A(2v^), (A.IO) 
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where Iy{-) is the modihed Bessel function of the hrst kind. Closed-form expression of the 
inverse Laplace transform C~^[f] exists if we can ensure that 




9f(Q!-|-l-|—> 0 , 


(A.ll) 


where 9f(-) means taking the real part. It will be shown later that this requires choosing the 
plus sign for a in flA.Sp so that 

Letting p = At + Ct^ and applying the inverse transform to /, we have 

4I+C+77 l-T+ioo „/(p-Ai) u(p-At) 

= / e ^tP dp 

^iTTlCyt Jt—100 

A+c+fi 


27riC'f 


3_(2-|- ^ 

Af u+Aju / AtU 


a 


-e 


e^tp ^ dp 

2 




Expressed in terms of v and v\ g is fonnd to be 


= e‘-C'-«^exp |' ' 


This completes the proof. 


An alternative method for solving (I A. 21) 


Instead of performing the Laplace transform of / with respect to u' and taking advantage of 
the Riccati system of ODEs, one may solve flA.2p nsing the method of characteristics. To 


proceed, we let / = / exp ( (a — 1) g 6^5 ds j, snch that / satishes the following PDE: 


_ d{{etu + Af) , ^ ^ Q 

2 du^ du dt 

where (3 = —e'^a — k. The terminal condition is given by 


(A.12) 


fit', u; t', CO, p, u') = 5{u — u). 


By applying the Laplace transform to / with respect to u\ 


Jo 
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the second order PDE flA.12j) is reduced to a first order linear PDE: 


f - - 9.5) f - «/ = 0. 


(A. 13) 


subject to 


Lemma 2. The solution to /iA.13\) is given by 

u'A,i \ ( 1 


f{t,^]t',u,r],u') = exp 


1 + Ct^ / \ 1 + Ct^ 


2/3/e" 


(A. 14) 


where 




^ f ds. 

2 Jt 


Proof. Let C, = ^At and f(tX',t',u;,ri,u') = f(t,^;t',u;,ri,u'). Then, flA.lSD can be written as 


dt 2AtdC ^ a /’ 

with terminal condition: /(t', C; t', cu, rj, u') = e~'^'^. 


(A.15) 


Notice that flA.lSp is a hrst-order PDE and its characteristic equations are given by 


d/ dC dt 

AlVC/ -ATh\^/2 1 • 

If we dehne (fit, (, /) = In / + ^ In /j' and 'fit, Ci /) = ^ + y // ds, then according to the 
characteristic equations, we have 


d0 

dfj 


f C 


dC 


“ C? ■ = 


0 . 


As a result, the solution to flA.lSp must have the following form: 


F{(f,tf) = F 


In / + ^ In C, 


1 £2 
2 



0 , 


where F is a function to be determined by the initial condition. Equivalently, we could rewrite 
the general solution in an explicit form of /: 


In/ 


-^InC + F 
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where if is a function to be determined by the initial condition. 
Using the initial condition; In / = —n'C) obtain 


As a result, H is given by 




HI I “' 2/3 

H iz) = -In z. 


Therefore, the solution to flA.lSp is given by 

u'C 


-2 rt' 


- 2 / 3/£2 


f{tX]t',u,ri,u') = exp 


1 H- ( / A ^ ds 


By noting that ( = AX and Ct = yAf A^ ^ ds = y jj' obtain (IA.14p . □ 


Finally, / is obtained by taking inverse Laplace transform C ^[f] which can be expressed 
explicitly, provided that 

m{/3) > 0. (A.16) 

In this case, it requires choosing the minus sign for a in flA.Sp . so that 


a = -|^ + ^|-c. 


(A.17) 


Using formula (lA.lOp and letting p = 1 + CX, we have 
f{t,u-,t’,uj,r],u’) = 


f-T+ioo p-2c-l 

27ri .U_4„ Ct 


u{p- 1) - Atu'{l - 1/p)' 
exp ( - - - ) dp 


exp ( I ^ ^r+ioo ^ 

e^Ap dp 


a 


exp 


t 

u-\-Atu' 

Ct 


27ri 


Ct 




Expressed in terms of the original variables v and v', g is shown to have the same expression 
as fl3.6l) . 


The choice for value of a 


When dealing with the inverse Laplace transform of / given by either flA.QD or flA.ldp . we need 
some technical condition which imposes constraints on the choice for value of a. Here, we 
briefly discuss how different technical conditions lead to different choices for value of a. 





























Denote c = cr + 'ici. Since 





io; + 


2ir] 





we conclude that cr > 0 under the usual branching line convention. Moreover, 


c 


2 

R 







implying that 




In conclusion, we must have cr > ^ + -p. Consequently, we have to choose « = — (| + ^) + c 
in order to meet flA.llD and = + —cin order to satisfy flA.lbp . 


Appendix B The 3/2 plus jumps model 


In this section, we extend our results to the 3/2 plus jumps model, which is given by 


1 c 

-^ = (r-q-X^)dt+ ^/Vt{pdWl + Vl-p"dID/) + (e^ - l) dNt, 

hi 

dVt = VM - i^Vt) dt + dWl 


(B.l) 


where a jump component modeled by a compound Poisson process is appended to the asset 
price process of the original 3/2 model fl2.1l) . Here, Nt is an independent Poisson process with 
constant intensity A and J stands for the random jump size which is assumed to be a normal 
random variable with mean p and variance . Moreover, the additional term in the drift X'ddt 
is the compensator for the jump component, where 


•d = E - 1] = - 1. 


With jumps in the asset price, the quadratic variation process is now given by A = f* W ds + 
and the triple transition density G satishes the following PIDE: 


'W 


v\ dG vd^G dG 


dG 


g - A'd - - —-h -vrw + v— + v{6t - nv)— + 




d^G 


2/ dx 2 dx"^ 


dy 


dv 


2 dv^ 


+ pev^ 


d^G 

dxdv 


+ X (G(t, x + z,y + z‘^, V] t', x', y', v') - G{t, x, y, v, t', x', y\ v'))fj{z) dz, 


where fj = is the probability density of the jump size J. The partial transform 

dehned by fl3.2l) admits a similar form given by 


G{t, X, y, V] t\ w, p, v') = x; t', w, p, ■ 
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where g satisfies that following PDE: 


dt 


ico [r — q 




2 ^ 

■ 00 - + igv + exp 






2(1 - 2ir7CT2) ) 


+ [v{0t - 


2^^9 


Kv) + iupev ] — + 
ov 


^2^3 Q 2 g 

2 dv'^ 


A 

— 2ir]a‘^ 



9 


The above PDE can be solved in the exactly same manner and the solution has the same 
expression as in fl3.6p with a modified parameter a given by 


a = ia;(r — q — A'd) + exp 


2ip{oj + rjfi) — oo^a"^ 


2(1 - 2iga'^) J ^1 - 2ir]a^ 


-A. 


Appendix C Proof of Theorem [2] 


To determine the new measure under which the calculation can be done conveniently, we need 
the following lemma. Let D be an open subset and Yf G VUdV be a one-dimensional stationary 
time-inhomogeneous diffusion process satisfying 

dYt = p{Yt,t) dt + a{Yt)dWt, 


where Wt is a Brownian motion under the measure Q. Its infinitesimal generator is given by 

= piy,t)dy + ^a^iy)dl 

Lemma 3. Assume 7 ( 2 /) is twice differentiable and 7 ( 1 /) > 0 for any y E T>. Define a new 
measure Q by the Radon-Nikodym derivative 


dQ 


Tt 


'liYfA _ A c 9 i(Ys) 
—L—d-e h 7(ys) 


(C.l) 


Suppose Yt remains to be finite and the boundary dV is unattainable under both Q and Q, then 
Q is equivalent to Q and 





~iYs) 



l{Yt)E 




(C.2) 


for any bounded measurable function l{y), where ]E([-] is the expectation taken under Q. More¬ 
over, under the new measure Q, the dynamics ofYt is modified to be 


dY, 


piYt,t) + a^{Yt) 


liXt)) 


dt + aiY,) dm, 


(C.3) 


where Wt is a Q-Brownian motion. 
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Lemma [3] is an extension of a related theorem established by Hurd and Kuznetsov (2008) 
from a general time-homogeneous diffusion process to a time-inhomogeneous diffusion process. 
It can be proved in a similar way based on the integrated Ito formula and the application of 
the Girsanov theorem. Equipped with Lemma [21 we are now able to obtain the conditional 
characteristic function of the integrated variance in closed form. 


Note that V)] = and Et[e-^* Us can be evaluated in 

two different ways for any zu, where Ut is a CIR process specified by (13.111) . On one hand, we 
consider the calculation of Et[e-^‘ based on Lemma [3l We hrst need to hnd a twice 

differentiable function 'y{y) such that 


7{Us) Us' 


(C.4) 


for some constants C. It can be easily verihed that 7 ( 2 /) = y'^ is a solution. In fact, since 
y,{y, t) = {k, + e^) — 6ty and a{y) = —ey, we have 


U^l{Us) a , + + 

+- u. - 

Apparently, flC.4p holds if G = —OgV and v is the solution to 

^^^ + («^+y)^ + ie = 0. 


1 ) 


which admits 



Applying Lemma [3] to Ut with 7 ( 2 /) = y'^ and l{y) = e yields 




r 

gjt 


177 




cPjiUs) 

7(G) 


ds — 
e ^ 


^-u sf Os 

t' 


(C.5) 


where pu is the transition density of the CIR process Ut under the new measure Q. Moreover, 
the Q-dynamics of Ut is given by 


dUt = {{k + e" + eV) - OtUt) df - e^/WtUWl, 


(C.6) 


where W} is a Q-Brownian motion. According to Lemma [H the transition density of Ut under 
Q is given by 


PuiUt'\Ut) = ^exp 


AtUt> + Ut\ fAtUt' 


Ut 


Ct 




A+%+ 2 i/ 
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. (C.7) 


















Notice that in order for make zero an unattainable boundary of Ut under Q, the Feller condition 
+ e^v) > £^/2 must be observed. Then, it follows that 


K 1 

£2 2 


1 

^ o 1 o • 


2 ' 


(C.8) 


On the other hand, we have the following alternative representation: 


E, 


e-h* 


= E* 






// 7^ 

gd t Us 


f; d. 

gJ t Us 


Uf 

Ut' 


By equating flC.SD and fIC.QF we obtain 


poo 

if TT- ds 

Us 

_ 1 

Jo 




Pu{Ut'\Ut)dUt'. (C.9) 


Pu{Ut'\Ut)dUt'. 


Noticing that the above equation holds for any w, we can then take the inverse Fourier trans¬ 
form with respect to w on both sides to obtain 


E+ 


d. 


Ut' 


Aifc) Pumut) 


Pu{Ut'\Ut) 

Substituting fl3.12p . flC.7p and fIC.Sp into fIC.lOp . we have 


(C.IO) 


E. 


/; ds 

gJ t Us 


Ut' 


h+^ ^tUfUt'^ 


This immediately leads to fl3.14p . 


Appendix D Proof of Corollary [T] 

Apparently, we have 

poo poo 

E^[gi^-Y,,+ia;/,,] ^ y^. ^ ^ ^lUtXt+ivh / y^. 

Jo Jo 

Then, it suffices to compute 

poo 

h{t,v;t',uj,p) = / g{t,v;t',u,p,v')dv' 

Jo 


+' +-1 -At 


a 


t Jo 


_u+A^ !' AtU'\ ^ ^2 

a 


= d«' 


u 


-^ ^ ) dz, 

Jo \ Ct J 


/ 2y/uz \ 
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where = Atu'. Using the equation 


e-‘t‘I,(\Vi) dt = ,p, X), 


r(^i) s 


where (j)=l + i + q/2, ^jJ = l + (^, X = ^ and 9^(0, s) > 0, we obtain 


h{t,v;t',u,r]) = 


-t)-x r(i-a) f u 


r(2c + i) \Ct 


M 1 - a,2c+ 1, 


u 

a 




,_.r(/3-d) / u 


r( 0 ) 

(/3 - a) f 
r(/3) \CtV 


u 






C,v 


Here, 


\ Hj ~ 

a = -^ + c, and B = 1 + 2c. 

2 


Note that the second equality follows from the identity: 


M(a, b, z) = e^M{b — a, b, —z). 
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